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Abstract 

Among the theories beyond the Standard Model (SM) of particle physics 
Supersymmetry (SUSY) provides an excellent dark matter (DM) candi- 
date, the neutralino. One clear prediction of cosmology is the annihi- 
lation cross section of DM particles, assuming them to be a thermal 
relic from the early universe. In most of the parameter space of Su- 
persymmetry the annihilation cross section is too small compared with 
the prediction of cosmology. However, for large values of the tan /3 pa- 
rameter the annihilation through s-channel pseudoscalar Higgs exchange 
yields the correct relic density in practically the whole range of possible 
SUSY masses up to the few TeV range. The required values of tan f3 are 
typically around 50, i.e. of the order of the top and bottom mass ratio, 
which happens to be also the range allowing for Yukawa unification in a 
Grand Unified Theory with gauge coupling unification. 
For such large values of tan f3 the associated production of the heavier 
Higgses is enhanced by three orders of magnitude and might be observable 
as one of the first hints of new physics at the LHC. 

1 Introduction 

Cold Dark Matter (CDM) makes up 23% of the energy of the universe, as deduced from the tem- 
perature anisotropies in the Cosmic Microwave Background (CMB) in combination with data on the 
Hubble expansion and other observations [I]. One of the most popular CDM candidates is the neu- 
tralino, a stable neutral particle predicted by Supersymmetry (SUSY) [21 OH]. The lightest neutralinos 
are spin 1/2 Majorana particles, which can annihilate into pairs of Standard Model (SM) particles. 
They are usually the Lightest Supersymmetric Particles (LSPs), which are stable if R-parity is con- 
served. This multiplicative quantum number forbids decays of SUSY particles with negative R-parity 
to known particles with positive R-parity, while the Next-to-Lightest SUSY Particles (NLSPs) and 



heavier ones will decay to the LSP and SM particles. 

SUSY is the most popular candidate for a theory beyond the SM, since it is the only theory so 
far, which solves several problems simultaneously: it provides a splendid candidate for DM, it has no 
quadratic divergencies in the Higgs sector, it allows for unification of gauge and Yukawa couplings 
and it predicts the Higgs mechanism by radiative corrections [5]. Search for SUSY particles is one of 
the prime objectives of the new Large Hadron Collider (LHC), which just starting taking data with 
several detectors jB]. But even if SUSY would be discovered at the LHC, it does not guarantee that 
the DM in the universe would be made of neutralinos, since many other candidates exist, see e.g. [7]. 

To prove that the DM is really made of neutralinos one would need to correlate additional prop- 
erties, like e.g. the mass or the annihilation cross section. The latter is a particularly good property, 
if DM is thermal relic from the early universe, since it is directly correlated with the relic density in 
an almost model independent way from the expansion history of the universe[3l Hj. Several studies 
exist on how well one can determine the annihilation cross section from accelerator data, see e.g. 
[8] and references therein. Up to now only methods in the so-called bulk region, focus point region 
or co-annihilation regions have been considered, since here rather specific SUSY mass spectra exist. 
In the so-called funnel region, where the annihilation proceeds through s-channel exchange of the 
pseudoscalar Higgs A, this has not been studied, since it is usually considered as a narrow stripe in 
parameter space for which the neutralino mass is close to the A-resonance, i.e. m x ~ m/i/2. Of course, 
one can consider rriA to be a free parameter and tune the mass to get the correct relic density or for 
that matter tune the masses of other SUSY particles contributing to annihilation. In this case the 
relic density does not constrain. However, in specific SUSY breaking scenarios all SUSY masses are 
related to each other via renormalization group equations [5]. In this Letter we consider the popular 
mSUGRA model with unified breaking scales for the fermions and bosons, respectively. This scenario 
allows for unification of gauge and Yukawa couplings and radiative electroweak symmetry breaking by 
minimizing the Higgs potential. In this case the mass of the pseudoscalar Higgs is sensitive to tan/3, 
the ratio of vacuum expectation values of the neutral components of the two Higgs doublets. Since 
tan/3 is a free parameter one expects that the correct relic density or equivalently the annihilation 
cross section can be obtained in almost any region of parameter space by tuning tan/3. It is the 
purpose of the present letter to see if this is true and discuss this in the context of future LHC physics. 

As it turns out, in a large region of parameter space the value of tan/3 is required to be around 
50, which implies an enhancement of the pseudoscalar Higgs production at the LHC by three orders 
of magnitude, since this cross section is cx tan 2 /3[9]. It should be noted that for these large values of 
tan /3 the contribution from other annihilation channels can be neglected, so discovering the Higgs at 
the LHC allows for the possibility to determine the relic density from LHC data and thus establishing 
a connection between the DM in the universe and the neutralinos produced in laboratory experiments. 
The value of 50 is furthermore of interest, since it is of the order of the ratio of top- and bottom-quark 
masses, which occurs naturally in SO(10) GUT theories. They feature in addition gauge coupling 
unification [TU1 [TT], Yukawa coupling unification [J3] and right-handed neutrinos, which are needed 
to generate small neutrino masses via the see-saw mechanism [14]. 

Several heavy flavour physics observables are enhanced by large values of tan /3 as well (see e.g. [TS] 
and references therein). Especially, the FCNC decay of B s — > fjbfj,, which proceeds in the SM via loops 
involving top quarks and W-bosons, is enhanced in SUSY by tan 6 /3. Experimental constraints from 
this tan /3 sensitive variable will be discussed, but excluded regions in SUSY parameter space by other 
low energy observables are outside the scope of this paper, since deviations from the SM are at the 
2-3<r level. The errors are dominated by statistical and systematic uncertainties, so the excluded or 
preferred regions depend on the assumed errors and the applied statistical treatment, which results in 
large differences between the different analysis, see e.g. for some recent analyses Refs. |16 | 1171 \TE[ IT9] 
and references therein. As said, discussing these analyses is beyond the scope of this letter, which 
concentrates on the SUSY parameter space allowing for a DM relic density determination from Higgs 
production at the LHC. The precision with which this can be done depends on the LHC luminosity 
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Figure 1: Annihilation diagrams for the lightest neutralino, which is a linear combination of the 
gaugino and Higgsino states: \x >= N\\B > +N 2 \Wq > +N 3 \Hi > +N 4 \H 2 >. The 
dependence of the amplitudes on masses and neutralino mixing parameter iVj has been 
indicated. 



and will be studied in a future paper. 

2 Predictions from mSUGRA 

The relic density and annihilation cross section a are related through: 

nh 2 = (i) 

< av > 

where the annihilation cross section a averaged over the relative velocities of the neutralinos is given 
in pb [3] and h rj 0.71 is the Hubble constant in units of 100 (km/s)/Mpc^. The best value for the 
relic density is Qh 2 = 0.1131 ± 0.0034 pQ. For the relation between cross section and expansion rate 
in Eq. [I] one assumes DM was a thermal relic, which froze out at the time when the annihilation rate 
was about equal to the expansion rate, given by the Hubble constant. For a given relic density Q 
the annihilation cross section is known independent of a specific model, since it only depends on the 
observed Hubble constant and the observed relic density. Its value is furthermore largely independent 
of the neutralino mass m x (except for logarithmic corrections) [3j H] . The DM constraint should exist 
for any model, but to be specific the mSUGRA model, i.e. the Minimal Supersymmetric Standard 
Model (MSSM) with supergravity inspired breaking terms, will be considered [5]. It is characterized 
by only 5 parameters: m,Q, fn\i2i tan/3, sign(/i), Aq. Here vtlq and my^ are the common masses 
for the gauginos and scalars at the GUT scale, which is determined by the unification of the gauge 
couplings at this scale. Gauge unification is perfectly possible with the latest measured couplings at 
LEP |llj . We only consider the dominant trilinear couplings of the third generation of quarks and 
leptons and set its common value at the GUT scale Aq equal zero, but the low energy values are 
different for all generations because of the different radiative corrections. Aq values different from zero 
were found not to influence the relic density, but are important if the relic density is considered in 
combination with B s — > fi^. Electroweak symmetry breaking (EWSB) fixes the scale of \i [5], so only 
its sign is a free parameter. The positive sign is taken, as suggested by the small deviation of the SM 
prediction from the muon anomalous moment [12J. 

The Born-level neutralino annihilation diagrams are shown in Fig. [TJ The cross sections are 
proportional to the final state fermion mass, which originates either from the Yukawa couplings for 
the Higgs exchange diagram or from the helicity suppression at the low energies involved in cold DMA 
|20j. Therefore heavy fermion final states, i.e. third generation quarks and leptons, are expected to 
be dominant. The Higgs exchange diagram is in addition proportional to tan/3 for down type quarks 
and l/tan/3 for up type quarks, so top quark final states are suppressed for large tan /3 and the main 
branching ratios will be into b-quarks, tau-leptons and muons. The W- and Z-final states have usually 




Figure 2: Left:Running of Higgs mass terms for different values of tan/3 for niQ = 1450 and m^ 2 = 
175 GeV. One observes that for large tan /3 the pseudoscalar mass, which at tree level equals 
mf+m 2 ,, becomes small for large tan /3, so for large tan /3 the dominant annihilation channel 
is via pseudoscalar Higgs exchange. Right: The strong dependence of the relic density as 
function of tan j3. 



a much smaller cross section due to the weak couplings involved, as was explicitly calculated with the 
CalcHEP program[2T]. In the MSSM the two mass terms of the two Higgs doublets, denoted by m<i 
and mi, receive radiative corrections from the top- and bottom Yukawa couplings, respectively [5]. 
The common value of these masses at the unification scale is \J m 2 , + (i^, but the large top Yukawa 
coupling can drive m 2 . negative at the electroweak scale for heavy top masses, thus leading to radiative 
electroweak symmetry breaking (EWSB). The top mass was predicted in this way to be heavy well 
before the top discovery [22] . For values of tan/3 ~ mt(mt)/mb(mt) ~ 160/3 ~ 50 the bottom Yukawa 
coupling becomes of the same order as the top Yukawa coupling, which could lead to unification of 
all Yukawa couplings of the third generation[13j. But in this case also mi is driven close to zero as 
well or becomes negative, as demonstrated on the left-hand side of Fig. [2} This implies that the 
pseudoscalar Higgs mass, which at tree level is given by m\ = m\ + m 2 , becomes light for large tan/3. 
Given that the relic density by A exchange is proportional to l/(4m 2 — m A ) 2 and ttia decreases fast 
with increasing tan/3 the relic density is a strong function of tan/3, as demonstrated on the right-hand 
side of Fig. [2} The value of tan (3 to obtain the WMAP relic density was calculated for each point in 
the rriQ — m^ 2 plane by simply scanning tan /3. The relic density was calculated with the public code 
micrOMEGAs 2.4|23j using Suspect 2.41 as mSUGRA mass spectrum calculator [24J. 

Figs. [3] and |4] shows that tan/3 can be adjusted in the whole mo — plane to get a correct 
relic density: the right-handed panels show the relic density and the values of tan /3 needed for the 
correct relic density. As mentioned in the introduction low energy SUSY masses are excluded, since 
else the effect of SUSY in loops would have shown up at LEP or heavy flavour observables. Since 
Br(B s —J- fifi) is proportional to tan 6 /3 it is an important constraint for large tan/3 scenarios. But 
this observable is also sensitive to the trilinear coupling at the GUT scale Aq and if Aq is left free, 
both the relic density constraint and the 90% C.L. upper limit of Br(B s — > fifj,) < 4.7 • 10~ 8 can 
be met in the whole plane and no excluded region from Br{B s — > fi/j,) is observed. 

The values of uia for the correct relic density are shown on the left-hand panel of Fig. |4j One 
observes that for larger neutralino masses, i.e. larger values of m^, the values of uia increase, as 
expected, since the neutralinos have to be on the tail of the pseudoscalar Higgs resonance for a correct 
relic density, typically myi/2m x ~ 1.2 — 2.4. It is interesting to note, that the needed values of tan/3 
are in practically the whole region close to 50, as expected since in that case the Yukawa coupling of 
the bottom quark becomes of the same order as the one of the top quark, so both Higgs mass terms 
in Fig. [2] become small. 

In Fig. [3] the top left is excluded, since here the LSP is not a neutral particle, but the stau is 
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Figure 3: Left: The region for correct relic density for various values of tan /3. One observes that the 
"funnel" region varies with tan/3. The trilinear couplings are at the GUT scale and the 
sign of /i=+l. Right: Relic density in the mo — m^ plane after tuning the value of tan /3. 



the LSP. The bottom right is excluded, since here EWSB does not work, since u 2 becomes negative, 
where a is the mixing parameter in the Higgs potential. In the transition to these forbidden regions 
co-annihilation of the LSP with an almost mass-degenerate other SUSY particle is possible. The co- 
annihilation and other properties of these transition regions have been reviewed nicely in Ref. [26J and 
the reader is referred to this reference for details and original references. Here only the salient features 
of these regions important for this analysis are mentioned. In the transition to the forbidden region at 
the top left the staus and lightest neutralinos are nearly degenerate in mass, so they will freeze out in 
the early universe at the same temperature. This leads to co- annihilation of the stau and neutralino 
into tau leptons, thus reducing the relic density according to Eq. [TJ In the bottom right-hand corner 
a similar feature is observed: here EWSB does not work, while /i 2 becomes negative and the lightest 
chargino and neutralino become almost degenerate. This can be understood as follows: the starting 
point of the Higgs mass terms at the GUT scale (see Fig. [ij is \J m| + /i 2 , so for large values of mo, 
(i has to become small. Else the Higgs mass terms in Fig. [2] will not become negative before the 
electroweak scale, as required for EWSB. If fi is small, the LSP and chargino mass terms are both 
determined by the smallest diagonal matrix element in the mass matrices, which is (i in this case. 
So in the transition region chargino-neutralino become mass-degenerate and their co-annihilation into 
W-bosons dominates. Then tan/3 has to be reduced again in order not to have a too large annihilation 
cross section, as is clear from the low tan /3 regions in the right-hand panel of Fig. |4j 

Finally, at small values of mo and m^ 2 the sfermions are light and the left-hand diagram of Fig. 
[T] dominates and again tan (3 has to be reduced in order not to have a too large annihilation cross 
section. This region is traditionally called the bulk region, since it does not correspond to a narrow 
stripe. 

In summary, if one allows tan/3 to vary in the mo — mj^ plane, one obtains the observed relic 
density for any combination of mo and my2, he. the relic density allows all masses for the SUSY 
sparticles. Such an mSUGRA scenario with large tan /3 can be uniquely tested at the LHC because of 
the large cross section for heavy Higgs production, as will be discussed in the next section. 

3 Expected Constraints from Higgs Production at the LHC 

The MSSM has five Higgs particles: the light Higgs h with a mass below 130 GeV, two heavy neutral 
Higgses (a scalar and pseudoscalar one, denoted by H and A, respectively) and two charged Higgses[27J. 
If the heavy Higgses are well above the light Higgs, all masses of the other charged and neutral Higgses 
are nearly degenerate and the production cross sections for the heavy scalar and pseudoscalar Higgs 
are nearly identical, so in practice one obtains twice the cross section, since the interference between 
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Figure 4: The values of itia (left) and tan/3 needed for the correct relic density in Fig. 3 (right). 




Figure 5: Left: The pseudoscalar Higgs production cross section as function of tan /3, both for the 
gluon fusion diagram and associated Higgs production with a b-quark for the different Higgs 
mass values indicated. Since the b-quark production is mostly in the forward direction, 
the scale on the right-hand side indicates if at least one b-quark is required to be in 
the acceptance, defined by rj < 2.5, and have a transverse momentum above 20 GeV/c. 
Right: expected discovery reach for the ATLAS detector at 14 TeV and a luminosity of 
30 The region already excluded at the Tevatron [HU] and the expected exclusion 

reach after the initial 7 TeV run at the LHC[29j have been indicated as well (assuming a 
luminosity of 1 fb^ 1 ). These sensitivity projections for future LHC running of the ATLAS 
and CMS detectors are preliminary. 



5 



the diagrams with pseudoscalar and scalar Higgses is negligible, even if the decay is into the same 
final state. 

The cross section for Higgs production at the LHC has been extensively studied, both in the SM 
and MSSM. For references we refer to the excellent reviews by Djouadi[9J. For the SM the gluon fusion 
via a top loop is the dominant diagram. The next dominant diagrams are the annihilation of heavy 
quarks into Higgs bosons. In next-to-leading order these heavy quarks, mainly bottom quarks, are 
obtained from gluon splitting. However, in the MSSM the Higgs coupling to down(up)-type quarks is 
proportional to the tan ft{\/ tan ft), so the coupling to top quarks is strongly suppressed at large tan ft 
and the production via the top loop is not dominant anymore, but the production in association with 
b-quarks dominates. This tan ft dependence is shown on the left-hand side of Fig. [5] for various Higgs 
masses. The associated production cross sections were calculated at Born level with CalcHEP|21j. 
while the gluon fusion via a loop was calculated with the formulae from the Higgs Hunter's Guide|27j 
after correcting a misprint in the squark contribution. Note that for small tan ft the SM behavior with 
the production via the top loop dominating is recovered. 

For values of tan ft = 50 the cross section for associated production of b-quarks and Higgses is 
enhanced by a factor 2500, so Higgs bosons can either be easily found or excluded. Although the 
Higgs decay into tau final state has only a 10% branching ratio, it is still the preferred search channel 
because of the strongly reduced background. At the LHC one expects the preliminary discovery reach 
for pseudoscalar Higgses to reach 800 GeV, as shown in Fig. [5] on the right-hand side for a recent 
analysis by the ATLAS Collaboration |28j . Here a luminosity of 30 fb^ 1 at the maximum LHC energy 
of 14 TeV was assumed, which corresponds to several years of running. Similar results were obtained 
earlier by the CMS Collaboration [29]. The region excluded at present by the combined data from 
the CDF and DO Collaborations at the Tevatron (using 1.8-2.2 fb" 1 at a center-of-mass energy of 
1.96 TeV) has been indicated as well|30j. The expected exclusion for the initial LHC run at 7 TeV is 
indicated from the preliminary CMS study assuming a luminosity of 1 f 6 -1 [3Tj . 

If one compares the right-hand side of Fig. [5] with the rriA values plotted in Fig. [4] one sees that 
the Higgs search is sensitive to mi/ 2 scales up to 800 GeV. In mSUGRA the gluinos are 2.7m 1 / 2 > so in 
this scenario with the neutralino providing all the dark matter, the Higgs search is sensitive to gluino 
mass scales up to 2 TeV. 

4 Conclusion 

In this Letter it has been shown that if the neutralino makes up the dark matter of the universe, the 
neutralino annihilation cross section, as expected from cosmology for a thermal relic, requires a value 
of tan ft of the order of 50, i.e. of the order of the ratio of top and bottom quark mass in a large region 
of parameter space. 

For such large values of tan ft the associated production of the heavier Higgses is enhanced by 
tan 2 ft, i.e. by three orders of magnitude. Combining the expected reach of Higgs searches of about 
800 GeV with the cosmological preferred region, if one assumes the neutralino to be the dark matter 
candidate, leads to a discovery reach in the mSUGRA mass plane corresponding to gluinos up to 2 
TeV. Here we discussed the search for the pseudoscalar Higgs, but the heavy scalar Higgs has typically 
the same mass, thus enhancing the cross sections by a factor two, since the interference between the 
pseudoscalar and scalar Higgs production is negligible (zero at Born level). 

If only a fraction of the dark matter consists of neutralinos, the annihilation cross section increases 
according to Eq. [T] and so for a given neutralino mass the via mass has to be lower. So if via has 
to be lower for a given value of m\/2 one obtains an even larger reach in the viq — m^ 2 plane by the 
Higgs searches. 

It should be noted that we analyzed the DM constraint in the framework of mSUGRA. However, 
the DM annihilation depends mainly on the masses of neutralinos, the neutralino mixing, the pseu- 
doscalar Higgs mass and tan ft, if one is outside the bulk and co-annihilation regions. So in similar 
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scenarios with large tan /3, as expected from Yukawa unification, the relic density constraint can be 
fulfilled, which implies that the Higgs production could be the first sign of new physics beyond the 
SM. 
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